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Abstract 

The principle of detailed balance states that in equilibrium each elementary process is equilibrated by its 
reverse process. For many real physico-chemical complex systems (e.g. homogeneous combustion, het- 
erogeneous catalytic oxidation, most enzyme reactions etc), detailed mechanisms include both reversible 
and irreversible reactions. In this case, the principle of detailed balance cannot be applied directly. We 
represent irreversible reactions as limits of reversible steps and obtain the principle of detailed balance 
for complex mechanisms with some irreversible elementary processes. We prove two consequences of the 
detailed balance for these mechanisms: the structural condition and the algebraic condition that form 
together the extended form of detailed balance. The algebraic condition is the principle of detailed bal- 
ance for the reversible part. The structural condition is: the convex hull of the stoichiometric vectors of 
the irreversible reactions has empty intersection with the linear span of the stoichiometric vectors of the 
reversible reactions. Physically, this means that the irreversible reactions cannot be included in oriented 
cyclic pathways. 

The systems with the extended form of detailed balance are also the limits of the reversible systems with 
detailed balance when some of the equilibrium concentrations (or activities) tend to zero. Surprisingly, the 
structure of the limit reaction mechanism crucially depends on the relative speeds of this tendency to zero. 

Keywords: reaction network, detailed balance, microreversibility, pathway, irreversibility, kinetics 
PACS: 82.40.Qt, 82.20.-w, 82.60.Hc, 87.15.R- 



1. Introduction 

1.1. Detailed Balance for Systems with Irreversible Reactions: the Grin of the Vanishing Cat 

The principle of detailed balance was explicitly introduced and effectively used for collisions bv lBoltzmann 



(|1964| ). In 1872, he proved his i7-theorem using this principle. In its general form, this principle is formu- 
lated for kinetic systems which are decomposed into elementary processes (collisions, or steps, or elementary 
reactions). At equilibrium, each elementary process should be equilibrated by its reverse process. The ar- 
guments in favor of this property are founded upon microscopic reversibility. The microscopic "reversing of 
time" turns at the kinetic level into the "reversing of arrows" : the elementary processes transform into their 
reverse processes. For example, the reaction J2i otiM — > J2j ftjBj transforms into ^2j(3jBj — > £\ otiAi 
and conversely. The equilibrium ensemble should be invariant with respect to this transformation because 
of microreversibility and the uniqueness of thermodynamic equilibrium. This leads us immediately to the 
concept of detailed balance: each process is equilibrated by its reverse process. 

For a given equilibrium, the principle of detailed balance results in a system of linear conditions on 
kinetic constants (or collision kernels). On the contrary, if we postulate just the existence of an a priori 
unknown equilibrium state with the detailed balance property the n a system o f nonl inear conditions on 
kinetic constan ts appear. These conditions were introduced in by IWegscheider (1901) and used later by 



Onsager (|l93lh . They are known now as the Wegscheider conditions. 



For linear kinetics, the Wegscheider conditions have a very simple and transparent form: for each 
oriented cycle of elementary processes the product of kinetic constants is equal to the product of kinetic 
constants of the reverse processes. 



Email address: agl5301e.ac.uk (A. N. Gorban) 



Preprint submitted to Chemical Engineering Science 



August 10, 2011 



However, many mechanisms of complex chemical and biochemical reactions, in particular mechanisms 
of combustion and enzyme reaction, include some irreversible (unidirectional) reactions. In many cases, 
complex mechanisms consist of some reversible and some irreversible reactions, equilibrium concentrations 
and rates of reactions become zeroes, and the standard forms of the detailed balance do not have a sense. 

In physical chemistry, the feasibility of a reaction depends on the energies and entropies of system states, 
initial, final, and transition ones. Nevertheless, some combinations of irreversible reactions are impossible 
irrespective of the values of thermodynamic functions. Since Wegschcider's time it is known that the cyclic 
sequence of irreversible reactions (the completely irreversible cycle) is impossible. It is forbidden by the 
principle of detailed balance. In a similar way, the reaction mechanism A ^ B, A — > C , C — > B is forbidden 
as well as A ^ B, A ^ C, C ->■ B. 

Two fundamental problems can be posed: 

(1) Which mechanisms with irreversible steps are allowed, and which such mechanisms are forbidden 
by the principle of detailed balance? 

In accordance with our knowledge, this question was not answered rigorously and the general problem 
was not solved. Beside that, the procedure of determining the forbidden mechanisms was not described. 

(2) Let a mechanism with some irreversible steps be not forbidden. Do we still have some relationships 
between kinetic constants of this mechanism? 

In our paper, both problems are analyzed based on the same procedure. Substituting the zero kinetic 
constants by small, however not zero values we return to the fully 'reversible case', in which all steps of 
the reaction mechanism are reversible. Then, we analyze a limit case, in which small kinetic parameters 
tend to reach 0. 



Such an idea was applied previously to several examples. In particular, IChul (|197lh used this idea for a 
three-step mechanism, demonstrating that the mechanism A^B,A~^C,B~^C can appear as a limit 
of reversible mechanisms which obey the principle of detailed balance, whereas the system A ^ B, A — > C, 
C —± B cannot appear in such a limit. However, this approach was not applied to the general analysis of 
multi-step mechanisms, only to a few systems of low dimensions. 

Since Lewis Carroll's "Alice's Adventures in Wonderland" , the Cheshire Cat is well known, in particular 
its inscrutable grin. Finally this cat disappears gradually until nothing is left but its grin. Alice makes a 
remark she has often seen a cat without a grin but never a grin without a cat. 

The detailed balance for systems with irreversible reactions can be compared with this grin of the 
Cheshire cat: the whole cat (the reversible system with detailed balance) vanishes but the grin persists. 

1.2. Detailed Balance: the Classical Relations 

First, let us consider linear systems and write the general first order kinetic equations: 

Pi = ^{kjPj - kjiPi) . (1) 

3 

Here, pi is the probability of a state Aj [i = 1, . . . , n) (or, for monomolccular reactions, the concentration 
of a reagent A)- The kinetic constant kij > (i ^ j) is the intensity of the transitions Aj — > Ai (i.e., 
kij is k^j). The rate of the elementary process Ai — > Aj is kjiPi. The class of equations (fT]) includes the 
Kolmogorov equation for finite Markov chains, the Master equation in physical kinetics and the chemical 
kinetics equations for monomolccular reactions. 

Let > be a positive equilibrium distribution. According to the principle of detailed balance, 
the rate of the elementary process A; — > Aj at equilibrium coincides with the rate of the reverse process 
Ai 4 Aj : 

hjpf = %p- q . (2) 

For a given equilibrium, p° q , the principle of detailed balance is equivalent to this system ^ of linear 
equalities. To find the conditions of the existence of such a positive equilibrium that (|2|) holds, it is 
sufficient to write equations ([2]) in the logarithmic form, lnp° q — hipj q = hi kij — In kji, to consider this 
system as a system of linear equations with respect to the unknown lnp^ q , and to formulate the standard 
solvability condition. 

After some elementary transformation this condition gives: a positive equilibrium with detailed balance 
((2|) exists if and only if 

1. If k^ > then kji > (reversibility); 
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2. For each oriented cycle of elementary processes, Ai 1 — > Ai 2 — > ...Ai q — > the product of the 
kinetic constants is equal to the product of the kinetic constants of the reverse processes: 

i i 

1 I kj+iij = 1 I ( 3 ) 
3=1 3=1 

where the cyclic numeration is used, i q +i = i\. 
Of course, it i s suff icient to use in ^ a basis of independent cycles (see, for example the review of 



Schna kenbcrg ( 197 



Let us introduce the more general Wcgscheider conditions for nonlinear kinetics and the gene ralized 
mass action law. (For a more detailed exposition we refer to the textbook of Yablonskii et al ( 199ll) .) The 
elementary reactions are given by the stoichiometric equations 

► ^ hj-^j (r = l,...,m), (4) 

i 3 

where Ai are the components and a r i > 0, f3 r j > are the stoichiometric coefficients. The reverse reactions 
with positive constants are included in the list (01 separately. We need this separation of direct and reverse 
reactions to apply later the general formalism to the systems with some irreversible reactions. 

The stoichiometric matrix is T = (7 r i), j r i = firi ~ ctri (gain minus loss). The stoichiometric vector ^ r 
is the rth row of T with coordinates = f3 r i — a r i . 

According to the generalized mass action law, the reaction rate for an elementary reaction §Q is 

n 

W r = K Y\_ a i" ) ( 5 ) 

i=l 

where a{ > is the activity of Ai . 

The list (Q| includes reactions with the reaction rate constants k r > 0. For each r we define fc+ = k r , 
w+ = w r , k~ is the reaction rate constant for the reverse reaction if it is on the list ^ and if it is not, 
w~ is the reaction rate for the reverse reaction if it is on the list (UJ) and if it is not. For a reversible 
reaction, K r = fc+ /k~ 

The principle of detailed balance for the generalized mass action law is: For given values k r there exists 
a positive equilibrium a° q > with detailed balance, — w~ . This means that the system of linear 
equations 

jriXi = In fc+ — In fc~ = In K r (6) 

i 

is solvable (xi = lna° q ). The following classical result gives the necessary and sufficient conditions for 
the existence of the po sitive equilibrium a° q > with detailed balance (see, for example, the textbook of 



Yablonskii et all (|1991I) ). 

Proposition 1. Two conditions are sufficient and necessary for solvability of (0): 

1. If fc+ > then k~ > (reversibility); 

2. For any solution X = (A r ) of the system 



^i.e. ^rlri = for all i\ 



XT = [i.e. 2^ A r 7rt = for all i ) 
the Wcgscheider identity holds: 



in 



J(k+) X ^Y[( k r) Xr - W 



r—1 r—1 



Remark 1. It is sufficient to use in (0) a basis of solutions of the system X G {A 1 , • • • , A 9 }. 

Remark 2. The Wegscheider condition for the linear systems is a particular case of the general 
Wcgscheider identity (0). Therefore, the solutions X of equation ^ are generalizations of the (non- 
oriented) cycles in the reaction networks. The basis of solutions corresponds to the basic cycles. This basis 
is, obviously, not unique. 
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Remark 3. In equation unknown x% = lna^ are independent variables and vector x can take any value 
in R n . In practice, this is not always true. For example, for heterogeneous systems with solid components 
some activities may vary in a narrow interval or may be even constant (see the more detailed discussion 
below in Section ] 3. 5)) . We do not study multiphase equilibiria in our paper. 

Remark 4. All the closed chemical systems have linear conservation laws: conservation of mass, various 
sorts of atoms, electric charge and other conserved quantities. They are linear functions of the amounts Ni 
of chemical components Ai. There is a problem of uniqueness and existence of a positive equilibrium with 
detailed balance or without it for every set of values of the independent conservation laws. To solve this 
problem we need some properties of the connection between activities and concentrations, (ai) O (cj). We 
do not assume any hypothesis about this connection and study just existence of a positive equilibrium with 
detailed balance in the space of activities. The Wegscheider identity ft8ty gives a necessary and sufficient 
condition for this existence. 

In practice, very often k~ = for some r, whereas k+ > 0. In these cases, the standard forms of the 
detailed balance have no sense. Indeed, let us consider a linear reversible cycle with an irreversible buffer: 

A x ^ A 2 ^ . . . A n ^ Ai -> A) • 

This system converges to the state where only po > and pi = for i > 0. In this state, trivially, 
u>+ = w~ = and it seems that the standard principle of detailed balance does not imply any restriction 
on the kinetic constants. Of course, this impression is wrong. 

Let us consider this system as a limit of the system with a reversible buffer, A x ^ A n (both reaction rate 
constants are positive), when the constant of the reverse reaction is positive but tends to zero: fc 1< _ — ► 0, 
fci^o > 0. For each positive value fci<_o > the condition of detailed balance tu+ = w~ gives the 
Wegscheider identity ([3]) for the cycle A\ ^ A 2 ^ - --An ^ A\\ The product of direct reaction rate 
constants is equal to the product of the reverse reaction rate constants. This condition holds also in the 
limit fci^-o - ► 0. So, any practically negligible but positive value of the reverse kinetic constant implies the 
nontrivial Wegscheider condition for the other constants. 

If we assume that the negligible values of the constants should not affect the kinetic systems then 
this Wegscheider condition should hold for the system with fully irreversible steps as well. Therefore, the 
following way for the formalization of the principle of detailed balance for irreversible reactions is proposed. 
We return to reversible reactions in which the principle of detailed balance is assumed by the introduction 
of small k~ > 0. Then we go to the limit k~ —> (k~ > 0) for these reactions. 

It is worth mentioning that the free energy has no limit when some of the reaction equilibrium constants 
tend to zero. For example, for the ideal gases F = J2i^i {RT In + /U? — RT) , where Cj is the concentration, 
Ni is the amount and [i® is the standard chemical potential of the component Ai. In the irreversible limit 
some /i? — > co. On the contrary, the activities 

p(«g4) (9) 

remain finite (for the ideal gases, for example, a, = c;) and the approach based on the generalized mass 
action law and the equations = w~ can be applied in the irreversible limit. 

Below, we study systems with irreversible reactions as the limits of the systems with reversible reactions 
and detailed balance, when some reaction rate constants go to zero. We formulate the restrictions on the 
constants in this limit and find the finite number of conditions that is necessary and sufficient to check. 
First of all we have to discuss the necessary notion of cycles for general reaction networks. 



1.3. Main Results 

We develop three approaches to the detailed balance conditions for the systems with some irreversible 
reactions. The first and the most physical idea is to consider an irreversible reaction as a limit of a 
reversible reaction when the reaction rate constant for a reverse reaction tends to zero. The limits of 
systems of reversible reactions with detailed balance conditions cannot be arbitrary systems with some 
irreversible reactions and we study the structural and algebraic restrictions for these systems. 

The second approach is based on the technical idea to study the limits of the Wegscheider identities 
(JSj) . Here, it is very useful to apply the concept of t he general (no n linear ) irreversible cycles or path ways 
developed recently far enough for our purposes by ISchuster et all ( 200fjf) : iGagneur fc Klamt ( 2004h and 
other. Let us write all reactions separately (including direct and reverse reactions) ([4]). The general 
oriented cycle or pathway is a relation between vectors 7 r with non-negative coefficients : ^ r A r 7 r = 0, 
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A r > and ^ r A r > 0. For each system with all reversible reactions and detailed balance the Wegscheider 
identity ([SJ holds for any oriented cycle. Therefore, if an oriented cycle persists in the limit with some 
irreversible reactions, then, for A r > 0, the rth reaction should remain reversible and for this cycle the 
Wegscheider condition persists. 

This property motivates the definition of the extended (or weakened, Yablonskv et al ( 2010l )) form 
of detailed balance in Section 13.11 through the general oriented cycles and the Wegscheider conditions. 
Theorem Q] states that a system satisfies the extended form of detailed balance if and only if it is a limit 
of systems with all reversible reactions and detailed balance. One part of this theorem ("if) is proved 
immediately in Section 13.11 the proof of the second part ( "only if ) exploits the third approach and is 
postponed till Section [4] 

The third idea is to study the limits when some equilibrium concentrations (or, more general, activities) 
tend to zero. For systems with all reversible reactions, we can explicitly express the constants of the reverse 
reactions through the constants of the direct reactions and the equilibrium activities: just use the detailed 
balance conditions, w + (a cq ) = w~ (a cq ). Here, instead of 2m parameters, kf 1 (m is the number of reactions) 
we use m + n parameters: m reaction rate constants k^ and n equilibrium activities a^ q . In this description 
of the reversible reactions, the principle of detailed balance is trivially satisfied. Some reactions become 
irreversible in the limits when some of the equilibrium activities tends to zero. Surprisingly, the structure 
of the limit reaction mechanism crucially depends on the relative speeds of this tendency to zero. 

In Section SI we assume that a® q = consti x e 5i and study the limits e —> 0. The n-dimensional space 
of exponents S = (Si) is split by m hyperplanes (j r , S) = on convex cones. Each of these cones is given 
by a set of inequalities (7r,<5) =0 (r = 1,. . . , m). In every such a cone, the limit reaction mechanism for 
e — > is constant. 

Using this approach, we prove the second part of Theorem Q] and even more: if a system satisfies the 
extended form of detailed balance then it may be obtained in the limit e — > from a system with all 
reversible reactions with given fc+ and a° q = const; x e Si for some exponents <5,; (Theorem 2|) . So, all 
the three approaches to the consequences of the principle of detailed balance for the systems with some 
irreversible reactions are equivalent. 

The computational problem associated with the extended form of detailed balance is nontrivial. For 
example, the oriented cycles (pathways) form a convex polyhedral cone and we have to formulate the 
structural condition of the extended form of detailed balance for all extreme rays (extreme pathways) of 
this cone (Theorem [2]): if A r > for a vector A from an extreme ray then the rth reaction should remain 
reversible. Calcula t ion of all these ex t reme rays is a well known and computational expensive problem 
(|Fukuda fc Prodonl . Il996l : IPapin et al 120031: iGagneur fc Klamti 12004 . In Theorem d we significantly 
reduce the dimension of the problem. 

Instead of the set of all stoichiometric vectors -f r (r = 1, . . . , m) in the whole space of composition 
W 1 {n is the number of components, m is the number of reactions) it is sufficient to consider the set of 
the stoichiometric vectors of the irreversible reactions in the quotient space W 1 / S, where S is spanned by 
the stoichiometric vectors of all reversible reactions. The simple exclusion of the linear conservation laws 
provides additional dimensionality reduction. The application of reduction methods is demonstrated in the 
case study in Section [331 

In Section 13.31 we formulate the main results for the simple case of linear (monomolecular) systems. 
Sections 13.41 and 13.51 are devoted to examples of nonlinear systems. In Section 13.41 the simple examples 
with obvious lists of the extreme pathways are collected. In Section [3.51 we analyze the possible irreversible 
limits for a complex reaction of methane reforming with CO2. 



2. Cycles and Pathways in General Reaction Networks 



There exist several graph representations of general reaction networks (lYablonskii et all . Il991t iTemkin et al , 



19961 ) and each of them implies the correspondent notion of a cycle. For example, each input and output 
formal sum in the reaction mechanism Q can be considered as a vertex (a complex) and then a reaction 
with the positive rate constant is an oriented edge. This g raph of the tra nsformation of complexes is 
convenient for the analysis of the com plex balance condition (lFeinberd . [l972l) . 

The bipartite graphs of reactions ( Volpert fe Khudvaev . 19851) gives us another example: it includes 
two types of vertices: components (correspond to Ai) and reactions (correspond to elementary reactions 
from Q ) . There is an edge from the ith component to the sth reaction if a r i > and from the sth reaction 
to the ith component if /3 r j > 0. The correspondent stoichiometric coefficients are the weights of the edges. 
This graph is convenient for the analysis of the system stability, for calculati o n of Jacobians for the right 
hand sides of the kinetic equations and for analysis of their signs (llvanova . 1979; iMincheva fc Roussell. 
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20071 ). For nonlinear systems, these graphs do not give a simple representation of the detailed balance 
conditions. 

We need a special notion of a cycle which corresponds to the algebraic relations between reactions. 
Let us recall that we include direct and inverse reactions in the reaction mechanism Q separately. Each 
solution of ([7]) may be represented as follows: 

A i x (YriOtuAi -> Ej/V 1 ^ 
A 2 x (V,,,„,l, -V, V 1 ,. 

(10) 

A m X £ 



i 

Here, a, = Es A s a s i = E s A s /3 s j. Therefore, we need the following definition of a cycle. 
Definition 1. ^4n oriented cycle is a vector of coefficients A ^ with all \ > £/ia£ satisfies h!0\ 



Remark 5 . Cycl es in catalysis and, especially, in biochemistry are called pathways \Schuster et a\ , \200A : 
Pavin et a\ , 200 A ). An oriented pathway is an oriented cycle from Definition^ An extreme (oriented) 



pathway is a direction vector of an extreme ray of the cone A + . A solution of equation {^) (a non-oriented 
cycle) is a non- oriented pathway. 

Qualitatively these concepts have been introduced in the early 1940s by Horiuti who applied them 



to heterogeneous catalytic reactions flHoriutil . 119731) . Horiuti used them to eliminate intermediates of 



the complex catalytic reaction by adding the steps of the detailed mechanism first multiplied by special 
coefficients. As result of such procedure, the chemical equation with no intermediates is obtained. 

All oriented cycles form the cone A + (without the origin). Extreme ray of a convex cone is a face that 
is, at the same time, a ray. Each ray may be defined by a directional vector A that is an arbitrary nonzero 
vector from this ray. The cone A + is defined by a finite system of linear equations j7]) and inequalities 
A r > 0. Therefore, it has a finite set of extreme rays. 

For integer stoichiometric coefficients, a S i, (3 S i, any extreme ray is defined by an uniform linear systems 
of equations with integer coefficients supplemented by the conditions A^ > and 0. Therefore, we can 
always select a direction vector with the integer coefficients. For each extreme ray, there exists a unique 
direction vector with minimal integer coefficients. 

For monomolccular reaction networks, these cycles coincide with the oriented cycles in the graph of 
reactions (where vertices are reagents and edges are reactions). 

There exists an oriented cycle of the length two for each pair of mutually reverse reactions. For these 
cycles the Wegscheider identities © hold trivially, for any positive values of k^. 

Remark 6. The systems without oriented cycles (A + = {0}) have a simple dynamic behavior. First of all, 
for such a system the convex hull of the stoichiometric vectors does not include zero: ^ conv{7i, . . . , 7 m }. 
Therefore, there exists a linear functional I that separates from {71, . . . , J m }: l(j s ) > for all s = 
1, . . . ,m. This linear function 1(c) increases monotonically due to any kinetic equation 

dc \ - 

-dt=^ Wsls 

s 

with reaction rates w s > 0: dZ(c)/df > if at least one reaction rate w s > 0. 



3. Extended Form of Detailed Balance 

3.1. Definition 

A practically important reaction mechanism may include reversible and irreversible steps. However, 
some mechanisms with irreversible steps may be wrong because they cannot appear as the limits of reversible 
mechanisms with detailed balance. Therefore, the first question is about the mechanism structure: what 
is allowed? 

The second question is about the constants: let the mechanism not be forbidden. If it is the limit of 
a system with detailed balance then the reaction rate constants may be connected by additional algebraic 
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conditions like the Wcgschcider conditions ©. Wc should describe all the necessary conditions. In this 
Section we answer both questions and formulate both conditions, structural and algebraic. 

We have to study study the identities © in the limit when some k~ — > 0. First of all, let us consider 
reversible reactions: if fc+ > then k~ > 0. It is sufficient to use in ([SJ only A with nonnegative 
coordinates, A r > 0. Indeed, the direct and reverse reactions are included in the list ([4]) under different 
numbers. Assume that X r < in an identity ([8]) for some r. Let the reverse reaction for this r have number 
r' . Let us substitute (kf) Xr in the left hand side of ((HJ by (kf,)~ Xr and (k~) Xr in the right hand side by 
(k~,)~ Xr . The new condition is equivalent to the previous one. Let us iterate this operation for various r. 
In the finite number of steps all the powers A r > 0. 

Let us use notation A for the linear space of solutions of ((TJ) and A + for the cone of positive solutions 
A (A r > 0) of 0. 

For reversible reactions, we proved the following proposition. Let the reactions are reversible and the 
direct and reverse reactions are included in the list (j4} separately. 

Proposition 2. The Wegscheider identity f3p holds for all A G A if and only if it holds for all positive 
A G A+. 

Elementary linear algebra gives the following corollary for reversible reactions. 

Corollary 1. The solution of the system of linear equations for logarithms of equilibrium activities (0|) 
exists if and only if for any positive solution A (X r > 0) of the system XT = the condition holds. 

Let us study identity (j8]) for a positive A when some of k r — > 0. In this limit, for every A G A + 
Corollary [T] gives two conditions: 

Corollary 2. Let k s > 0, k s — >• fc' lm and the Wegscheider identity holds for k s . Then 

1. If X s > and kf — > for some s then for some q X q > and k~ — > 0; 

2. If for all positive components X s > the limit constants are positive, fc' lm± > 0, then the condition 
(§) holds for k l ™ ± . 

We can interpret the positive solutions of ([7]) as oriented cycles (linear or nonlinear) . The first condition 
means that if a cycle is cut by the limit fc+ — > in one direction then it should be also cut by a limit 
k~ — > in the opposite direction: the irreversible cycle is forbidden. This remark leads to the definition of 
the structural condition of the extended form of detailed balance. 

Definition 2. A system of reactions Op satisfies the structural condition of the extended form of detailed 
balance if for every A G A + the reaction which satisfy X s > are reversible: if X s > then kf 1 > 0. 

This condition means that all cycles should be reversible. The second condition means that for all cycles 
A G A + which persist in the system with irreversible reactions the correspondent Wegscheider condition 
([8]) holds. This is the algebraic condition of the extended form of detailed balance. Now, we are ready to 
formulate the definition of the extended form of detailed balance. 

Definition 3. The subsystem satisfies the extended form of detailed balance if both the structural and the 
algebraic condition hold for all A G A+. 

The following theorem gives the motivation to this definition. 

Theorem 1. A system with irreversible reactions is a limit of systems with reversible reactions and detailed 
balance if and only if it satisfies the extended form of detailed balance. 

Proof. Let us prove the direct statement: if a system is a limit of systems with reversible reactions and 
detailed balance then it satisfies the extended form of detailed balance. Indeed, let a system of reactions 
be a limit of systems with reversible reactions and detailed balance. This means that for each j = 1,2,... 
a set of reaction rate constants kfj > is given, kfj > satisfy the principle of detailed balance for all j 
and 

kf = lim kf , . 

Assume that the structural condition is violated: there exists such a A G A + that X s > for an irreversible 
reaction (kf > 0, k~ = 0). For all j = 1, 2, . . . the principle of detailed balance gives: 

n ( fc ^) A " - n ■ (ii) 

r, \ r >0 r,\ r >0 



7 



If A r > then fc+ > 0. Therefore, for these r, sufficiently large j and some e,S > S > k^r- > e > 0. The 
left hand side of (fTTj) is separated from zero. The right hand side of (TTTI) tends to zero because all factors 
are bounded and at least one of them tends to zero, k~ ■ — > 0. This contradiction proves the structural 
condition. To prove the algebraic condition, it is sufficient to notice that the Wegscheider identity for 
kfj > holds for all j, hence, it holds in the limit j — > oo. 

We will prove the converse statement (if a system satisfies the extended form of detailed balance 
then it is a limit of systems with reversible reactions and detailed balance) in Section |4j in the proof of 
Theorem 01 □ 

3.2. Criteria 

All A G A + participate in the definition of the extended form of detailed balance. Nevertheless, it is 
sufficient to use a finite subset of this cone. 

We can check directly that if for a set {X s } the structural and the algebraic conditions of the extended 
form of detailed balance hold then they hold for any conic combination of {A 5 }, A = J2 s as ^ S > fls — 0- 
Therefore, it is sufficient to check the conditions for the directional vectors of the extreme rays of A+. 

Let a reaction mechanism satisfy the extended principle of detailed balance. If we delete from this 
mechanisms any irreversible elementary reaction or any couple of mutually reverse elementary reactions, 
the resulting mechanism satisfies the extended principle of detailed balance as well. 

A cone is pointed if the origin is its extreme point or, which is the same, this cone does not include a 
whole straight line. The cone A + is pointed because it belongs to the positive orthant {A | A > 0}. 

It is a standard task of linear progr amming and computational convex geometry to find all the extreme 
rays o f the polyhedral pointed cone A + ( Bertsimas fc Tsitsiklisl . [l997 ; Motzkin et al . 19531 iFukuda fe Prodon . 



19961 ). Let the directional vectors of these extreme rays be {A s | s = 1, • ■ ■ , q}. 



Theorem 2. The system satisfies the extended form of detailed balance if and only if the structural and 
algebraic conditions hold for the directional vectors {X s \ s = 1, • ■ ■ ,q} of the extreme rays of the cone A + . 

Theorem [5] follows just from the definition of extreme rays and the Minkowski theorem which states 
that every pointed cone given by linear inequalities admits a unique representation as a conic hull of a 
finite set of extreme rays. 

This criterion can be simplified as well: it is necessary and sufficient to check the structural conditions 
for the extreme rays of A + and then the algebraic condition for a maximal linear independent subset of 
{A s | S = l,.-. ,q}. 

Corollary 3. The system satisfies the extended form of detailed balance if and only if the structural 
conditions hold for all directional vectors {X s \ s = 1, • • • , q} of the extreme rays of the cone A + and the 
algebraic conditions hold for a maximal linear independent subset of {X s \ s = 1, • • • , q}. 

If, for a given reaction mechanism, the set {A s | s = 1, • • • , q} of directional vectors of the extreme rays 
of A + is known, then it is easy to check, whether this mechanism satisfies the structural conditions of the 
extended form of detailed balance. It is sufficient to examine for each A* > 0, whether k~ > 0. 

After these conditions are examined, it is a simple task to extract the independent set of the Wegscheider 
identities that should be valid: just select a maximal linear independent subset from the set of A s and 
write the correspondent Wegscheider identities. 

It is convenient to use all the extreme pathways especially if we would like to study all the subsystems 
of the given system, which satisfy the extended form of detailed balance. On the other hand, it is compu- 
tation ally expensive to find the set {A s | s = 1, •■ • ,q} (see, for example, the paper bv iGagneur fc Klamtl 
The amount of computation could be significantly reduced because it is not necessary to use all 



the extreme pathways. 

Let us consider a reaction mechanism, which includes both reversible and irreversible reactions. For 
the reversible reactions, let us join the direct and reverse reactions. Let 71, . . . ,7r be the stoichiometric 
vectors of the reversible reactions and V\, . . . , v s be the stoichiometric vectors of the irreversible reactions. 
We use r r for the stoichiometric matrix of the reversible reactions and A r for the solutions of the equations 
AIV = 0. 

The linear subspace S — span{7i, . . . , 7,,} C K™ consists of all linear combinations of the stoichiometric 
vectors of the reversible reactions. Let us consider the quotient space M. n /S. We use notation Vj for the 
images of i/j in IR' 1 /^- 

The following theorem gives the criteria of the extended form of detailed balance, which are more 
efficient for computations. 

Theorem 3. The system satisfies the extended form of detailed balance if and only if 
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1. The convex hull of the stoichiometric vectors of irreversible reactions does not intersect S, i.e. 



$. conv{z/i,...,i/ s }; (12) 
2. The reversible reactions satisfy the principle of detailed balance. 

Proof. Let the condition 1 be violated, i.e. S convjFi, . . . , T7 S }. In this case, there exist such a nonnegative 
9i > that Y^"j=i @j = 1 an( i Sj=i Qj v j £ ^- This means that $^ =1 ®i v j + Sl=i = 0- We can 
transform the sum Y^=i ^i7i m a combination with positive coefficients if for any negative Xi we substitute 
7i by the stoichiometric vector of the reverse reaction, that is, —7,. As a result, we get the element of 
A_|_, a combination of the stoichiometric vectors with nonnegative coefficients, which is equal to zero and 
includes some stoichiometric vectors of the irreversible reactions with nonzero coefficients. Therefore, the 
structural condition of the extended form of detailed balance is violated. 

Let the structural condition be violated. Then there exist a combination X)j=i + ^ili = 

with 9j > and Ylj=i ®i > ®- Let us notice that 

EOj Xi 

j=i ^i=i °J i=i ^=1 °o 

and, therefore, £ conv{77 1; . . . ,V S }. The condition 1 is violated. 

We proved that the condition 1 is equivalent to the structural condition of the extended form of detailed 
balance. 

If the condition 1 holds then the condition 2 is, exactly, the algebraic condition of the extended form 
of detailed balance. □ 

Remark 7. The first condition of Theorem, ^ conv{Fi, . . . , v s }, is equivalent to the existence of such a 
linear functional I on W 1 that l{vj) > for all j = 1, . . . , s and l{"fj) = for all j = 1, . . . , r. 

3.3. Linear Systems 

The r esults of previous Sect ions for a linear system ^ have a geometrically clear form (see also the 
paper bv lYablonskv et al ( 2010( )). 



Proposition 3. The necessary and sufficient condition for the extended form of detailed balance is: In 
any cycle Ai 1 Ai 2 As — > Ai 1 with the strictly positive constants ki j+l i< > (here i q +i = i\) all 

the reactions are reversible (ki j i j+1 > 0) and the identity (0) holds. 

The states (reagents) A q and A r (q ^ r) are strongly connected if there exist oriented paths both from 
A q to A r and from A r to A q (each oriented edge corresponds to a reaction with the nonzero reaction rate 
constant). From Proposition [3] we get the following statement. 

Corollary 4. Let a linear system satisfy the extended form of detailed balance. Then all reactions in any 
directed path between strongly connected states are reversible. 

In brief, a linear system with the extended form of detailed balance can be described as follows: (i) the 
oriented cycles are reversible and satisfy the classical condition ([3]), (ii) the system consists of the reversible 
parts and the irreversible transitions between these parts and (iii) the system of these irreversible transitions 
is acyclic. 

For example, let us analyze subsystems of the simple cycle, A\ ^ A2 ^ A3 ^ A\. 



r T 



-10110-1 

1-1 0-1 1 
1-1 0-1 1 



(13) 



The cone of nonnegative solutions A + to the equation XT — has extreme rays with the following 
direction vectors: A 1 = (1, 1, 1, 0, 0, 0), A 2 = (0, 0, 0, 1, 1, 1), A 3 = (1, 0, 0, 1, 0, 0), A 4 = (0, 1, 0, 0, 1, 0), and 
A 5 = (0,0,1,0,0,1). Vectors A 3-5 give trivial identities © k+k~ = k~kf (i = 1,2,3) and vectors A 1 ' 2 
give the same identity: k^k^k^ = k^k^k^. 

If we delete one elementary reaction from the simple cycle (i.e. one column from T T (|13|0 then one of 
the nonnegative solutions A 1 ' 2 persists and, due to the extended detailed balance principle, all the reactions 
should be reversible. This means that the structural condition of extended detailed balance is not satisfied 
for the simple reversible cycle without one direct or reverse reaction. If two reactions are reversible then 
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the third should be reversible or completely vanish. If we delete one direct reaction (with number 1 , 2 or 
3) and one reverse reaction (with number 4, 5 or 6) then there remain no non-trivial solutions in A + and, 
therefore, no non-trivial relations between the constants persist after deletion of these two reactions. 

For the linear systems, the oriented cycles in the graph of reactions (where vertices are the components 
and edges are the reactions) give the positive solutions to the equation (JT)): for a linear oriented cycle C 
the sum of the stoichiometric vectors of its reactions is zero. Moreover, any positive solution of ([7|) is a 
convex combination of such cyclic solutions and, therefore, the directed vectors of the extreme rays of A + 
can be selected in this form. 



3.4- Simple Examples of Nonlinear Systems 

In this section, we present several elementary examples. For these examples, the sets of the extreme 
pathways are obvious. 



ki &2 

Let us examine a reaction mechanism with irreversible reactions A — 5- B and 2B — ^> 2A. 



(14) 



The cone A + is a ray with the directional vector A = (2, 1). Both A1.2 > 0, hence, both reactions should 
be reversible and the condition holds: (kf) 2 k 2 = (kx) 2 k 2 ■ 

Let us slightly modify this example: 2H — > YL 2 , H + H2 — > 3H. 



-2 2 
1 -1 



(15) 



The cone A+ is a ray with the directional vector A = (1, 1). Both Ai,2 > 0, hence, both reactions should 
be reversible and the condition holds: k^k 2 = k^k 2 . 

If we change the direction of one reaction in the previous example then the new irreversible systems 
satisfies the extended form of detailed balance: 2H — > H2, 3H — > H + H2. 



-2 -2 
1 1 



(16) 



The cone A + is trivial (it includes no rays, just the origin), hence, the structural condition holds. The 
algebraic condition trivially holds, because there is no reversible reaction. 

Let us add the forth reversible and nonlinear elementary reaction A\ + A2 ^ 2A 3 (with the constants 
kf) to a linear reversible cycle. We should add to T T (fl"3)) two new columns: 



(17) 



The extreme rays of A + include four rays that correspond to pairs of mutually reverse reactions (A 1 4 ), 
two rays that correspond to the linear cycle (A 5 ' 6 ) and six rays for three nonlinear cycles (A 7 " 12 ): (i) 



Ai + A 2 ->■ 2A 3 , A 3 
A 3 ->■ A 2 , A 2 -» A x : 



A 2 , A 3 



Ai; (ii) Ai + A 2 -> 2A 3 , A 3 ->• Ax, Ax -> A 2 and (iii) Ax + A 2 -> 2A 3 , 
X 6 = (0,0,0,0,1,1,1,0), 



(1,1,1,0,0,0,0,0), 
(0,0,1,1,0,1,0,0), 



(0,1,0,0,0,0,1,1), 



A 9 = (1,0,2,1,0,0,0,0), A 10 = (0,0,0,0,1,0,2,1) 



A 11 = (0,0,0,1,1,2,0,0), 



(1,2,0,0,0,0,0,1) 



We omit A 1 4 which do not produce nontrivial conditions. For the reversible reaction mechanism (when 
kf_4 > 0), there are two independent Wegscheider identities (|8]) that formalize the classical principle of 
dcij cliLgc! bcilcLiiCG! k-y k^k^ := k-y k^ k^ ciiid k^ k^ k^ — — k^ k^ k-^ . If some of the elementary reactions are 
irreversible then the direction vectors A 5-12 produce 8 conditions. For A 5,7 ' 9,11 these conditions are below. 



• (A 5 ) If kf 2 3 > then k l 2 3 > and kfk^k^ = k 1 k 2 k 3 ; 



• (A ) If k 3i , k 2 > then k 3 4 , k 2 > and kfk^k 2 



fc 3 /c 4 k 2 ; 



• (A 9 ) If k+ 3 4 > then fcf 3 4 > and kf(k^) 2 kt = k^(k 3 ) 2 k^: 
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• (A 11 ) If fc+, fc^ 2 > then k^,k+ 2 > and k+k^{k 2 ) 2 = k 4 ;k+(k+) 2 . 

To obtain the conditions for .x 6 ' 8 ' 10 ' 12 it is sufficient to change the superscripts + to ~ and inverse. These 
8 conditions represent the extended form of detailed balance for a given mechanism. To check, whether a 
subsystem of this mechanism satisfies the extended form of detailed balance, it is necessary and sufficient 
to check these conditions. 



3.5. Methane Reforming Processes: a Case Study 
3.5.1. The System 

Methane reforming with CO2 is a complex reaction network (IBensonl . il 98 II ). The main reactions in the 
methane reforming arc: 

1. C0 2 + H 2 ^ CO + H 2 (RWGS, Reverse water-gas shift); 

2. CH 4 + C0 2 ^ 2C0 + 2H 2 (Dry reforming); 

3. C0 2 + 4H 2 ^ CH 4 + 2H 2 (Methanation); 

4. CH 4 + H 2 ^ CO + 3H 2 (Steam reforming); 

5. CH4 ^ 2H 2 + C (Methane decomposition); 

6. 2C0 ^ C0 2 + C (Boudouard reaction); 

7. C + H 2 ^ CO + H 2 (Coal gasification). 

For the reagents, we use the notations A\ = CH 4 , A 2 = C0 2 , A 3 = CO, A4 = H 2 , A 5 = H 2 0, ^6 = C. 
Amount of Aj is Ni. There exist three independent linear conservation laws: be = Ni + N2 + N3 + Nq; 
6h = 4iVi + 2N4 + 2N$; bo = 2N2 + Ns + N5. The number of degrees of freedom in the closed system is 
three (six components minus three independent conservation laws). 

This example enriches our discussion because it deviates from the nice abstract scheme discussed above. 
First of all, the reactions 1-7 are not elementary steps. We consider them as overall reactions which have 
their own intrinsic and complicated reaction mechanism. This does not cause a serious problem because 
the generalized mass action law describes the equilibria of the complex overall reactions as well as the 
equilibria of the elementary ones. Therefore, we can apply the concept of the extended form of detailed 
balance and our theorems [TH3] to the process network 1-7 build from the complex reactions. Rigorously 
speaking, we deal not with the elementary reaction steps but with the main equilibria and may discuss, 
for example, not the "Boudouard reaction" but the "Boudouard equilibrium" . 

The second problem is the heterogeneity of the system: A\ , . . . , A§ are gases and Ag — C is solid. Some 
of the reactions go on the surface of the solid. 

If a multiphase system is ideal and the solid components are stoichiometric ones (i.e. they have a fixed 
composition) then the free energy has the form 

F= N l (RT\nc l + $ - RT) + ^ . (18) 

Ai — gas At — solid 

Here, the free energy of solid components differs from the free energy of gases by the absence of the term 
RTNlnc. This term corresponds to the ideal gas pressure PV = NRT. In our case, 

5 

F = Y J Ni{RT\nc i + pP i -RT) + N & nl. (19) 
1=1 

To define the activities, we follow (j9]). For the ideal gases ai = a and for the stoichiometric solids a; = 1. 

In section \l.2\ we studied homogeneous systems and considered Xi = In a; as independent unknowns in 
the detailed balance equations ©: 

^ JriXi = In K r (xi = In a° q ) . (20) 

i 

Therefore, for any solution of this system, the activities ai = expxi represented a positive equilibrium. 

In a heterogeneous system with the free energy (|18[) the activities for the solid components are constant, 
the correspondent Xi = 0. Let x = (x;) be a solution to equations (|2T)|) , Tx = K, where K is the vector of 
the equilibrium constants. The vector a = (a^), = h\Xi is a vector of equilibrium activities if and only if 
xi = for all the solid components Ai. Instead of analyzing the solvability of the detailed balance equations 
(|20[) we have to study its solvability under additional condition: Xi = for all the solid components Ai. 

Let us postpone the discussion of the extended principle of detailed balance in multiphase systems and 
consider the system "gaseous mixture + one stoichiometric solid" . Let A n be solid. 
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If there is the only solid component then the solvability conditions for the system (|20p and for this 
system with additional condition x n = coincide. Indeed, there exist a positive stoichiometric linear 
conservation law: 



E 



7r 



for all r and bi > for all i . 



For example, this may be conservation of mass or of the amount of atoms. Let b = (bi). For any solution 
of the detailed balance conditions (|2T)|) x = (xj), the vector 



b n 

is also a solution to (|20|) with the condition x' n = 0. 

So, for our example with seven equilibria 1-7 the conditions of the extended principle of detailed balance 
for the heterogeneous system with solid A§ =C are described by the theorems [TH3] and we can use the results 
of the preceding sections. 



3.5.2. The Classical Wegscheider Conditions 

To formulate the classical Wegscheider identities, we have to join the direct and inverse reactions and 



to find the basic solutions of the system of linear equations Ar = 0. 
example is: 





The stoichiometric matrix for this 



-1 
1 

-1 
1 




-1 
-1 
2 
2 





1 

-1 


-4 
2 




-1 

1 
3 
-1 






1 

-2 


1 






1 

1 

-1 
-1 



(21) 



The system of seven equations XT = is redundant. There are only three independent equations (one 
equation for every degree of freedom) . It is sufficient to take the components of stoichiometric vectors that 
correspond to the components A 2l A 4 , Aq. Other components satisfy the same linear relations as the 
selected ones. The reduced matrix r 1 is 



-1-1-10 1 
-1 2 -4 3 2 1 
1 1 -1 



(22) 



There are four independent solutions of the equations XT = (seven variables minus three independent 
equations). For example, we can take the following basis of solutions: (—1, 1, 0, —1, 0, 0, 0), (0, 0, 0, —1, 1, 0, 1), 
(1,0,0,0,0,1,1), (1,0,-1,-1,0,0,0). 

The correspondent Wegscheider identities arc: K 2 — K1K4, K 5 K 7 = K4, KiKqKj = 1, K\ = K3K4. 



3.5.3. Allowed and Forbidden Mechanisms 

In general, all the seven reactions can be considered as reversible but under various conditions some 
of them are almost irreversible. Let us study which combinations of irreversible reactions are possible in 
accordance with the extended form of detailed balance. 

For example, existence of the positive solution (0,1,0,0,0,1,1) G A guarantees that the irreversible 
system C0 2 + H 2 -4 CO + H 2 0, CH 4 + C0 2 -4 2CO + 2H 2 , C0 2 + 4H 2 -4 CH 4 + 2H 2 0, CH 4 + H 2 
CO + 3H 2 , CH 4 -> 2H 2 + C, 2CO -> C0 2 + C, C + H 2 -> CO + H 2 is forbidden by the extended 
form of detailed balance. This conclusion is also obvious from the correspondent Wegscheider condition 
K2KqKj = 1. Indeed, if all the k~ — > for bounded from below > e > then all Kj — > 00 and 
K2KQK7 — > 00. This contradicts the Wegscheider condition. 

The firs t reaction (RWGS, Reverse water-gas shift) is reversible in the wide interval of conditions 
(jMoel . Il962h . Let us first study all the reaction mechanisms with the reversible first reaction and the 
irreversible reactions 2-7. We find the combinations of the directions of the irreversible reactions that 
satisfy the extended form of detailed balance. As a criterion of the extended form of detailed balance we 
use Theorem [3] After that, we consider other reactions as the reversible ones (in addition to RWGS) and 
study the corresponding reaction mechanisms. 

The space S is a straight line with the directional vector 71 with coordinates (—1, —1, 0) in the coordinate 
system (N2, iV 4 , Nq) that corresponds to the components A2, A±, Aq. Let us represent the quotient space 
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Figure 1 : Images of the stoichiometric vectors of irreversible reactions Vj = ±7^ in R J / S for various combinations of directions 
of reactions l|24|l in coordinates N2 (abscissa), Nq. The configurations with £ conv{!72, . . . ,1^7} are outlined. Vectors i7 2 , V3 
and F4 coincide as well as vectors and F7. 



IR 3 //^ in the coordinate system (N2,N$) that corresponds to the components A2, A$. For this purpose, we 
have to eliminate the coordinate iV 4 using vector 71. As a result, we get the following vectors: 



72 
75 



>7 3 



.7 6 = 



■ 74 



,77 = 



(23) 



For example, to find 7 2 , we take 72 (the second column in ([22])) and exclude the coordinate N4 by adding 
271. The result is a vector 72 + 271. In coordinates (N2, Ne), this vector gives us 7 2 . 

The stoichiometric vectors of irreversible reactions arc +7j or — jj (j = 2,..., 7). Their images in 
the quotient space R 3 /S* are +7^ or — 7^. The extended form of detailed balance requires that the convex 
envelope of these vectors should not include zero. We have to arrange signs in ±jj to provide this property. 
First of all, we see immediately from (|23[) that the second and the forth reaction should have the same 
directions and the third reaction should have the opposite direction. The directions of the sixth and the 
seventh reactions should be opposite. Therefore, we have to analyze eight possible reaction mechanisms. 
Let us represent them by the directions of reactions: 



(a) 



(b) 



(c) 



(d) 



(e) 



(f) 



(g) 



00 



(24) 



Arrows here correspond to the directions of reactions. For example, the case (a) corresponds to the reaction 
mechanism 



1. C0 2 



H 2 



CO + H 2 0; 



CH 4 + C0 2 - 
C0 2 + 4H 2 <- 
CH 4 + H 2 - 
CH 4 -> 2H 2 - 
2CO -)• C0 2 



► 2CO 
CH 4 - 

y CO 4 
C; 

-C; 



h2H 2 ; 
2H 2 0; 
3H 2 ; 



7. C + H 2 <- CO + H 2 . 

Combinations (c) and (f ) contradict the condition 1 from Theorem [3] the origin belongs to the convex 
envelope of the vectors Vj of irreversible reactions (see Fig. [1]). Hence, only six combinations of directions 
of irreversible reactions satisfy the extended form of detailed balance (from 2 6 = 64 possible combinations 
of directions): (a), (b), (d), (e), (g) and (h). 

Let us extend the list of reversible reactions. If we assume that the second reaction (dry reforming), 
is reversible together with the first one (RWGS) then the third and the forth reactions should be also 
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reversible because 73 = 271 — 72 and 74 = 72 —71, hence, 73.4 £ span{7i,72}. According to the condition 

1 from Theorem [31 this contradicts to the extended form of detailed balance if the first and the second 
reactions are reversible and the third and the forth are not. 

Analogously, in addition to the reversible reaction RWGS, the Boudouard equilibrium 6 and coal gasi- 
fication 7 can be reversible only together because 77 = —71 — 76. 
We have to consider three possible sets of reversible reactions: 

1. 1, 2, 3 and 4; 

2. 1 and 5; 

3. 1, 6 and 7. 

For all three cases, dim S = 2 and dim(IR 3 /S') = 1. We will use for the quotient space the coordinate 
Nq which corresponds to Aq = C. 

In the first case, let us exclude the coordinate A4 from 7 5 6 7 (|2"5)l using vector 7 2 . We get one- 
dimensional vectors 

75 = !, 7e = 7 7 = -1 ■ 

To satisfy the extended form of detailed balance the directions of the fifth and the sixth reaction should 
coincide and the direction of the seventh reaction should be opposite: there are two possible combinations 
of arrows in irreversible reactions 5, 6 and 7 if reactions 1, 2, 3 and 4 are reversible: 5 — >, 6 — >, 7 <— and 
5 6 7 -h 

In the second case, let us exclude the coordinate A4 from 723467 p3)l using vector j 5 . We get 
one-dimensional vectors 

7 2 = -2/3, 73 = 2/3, 74 = -2/3, 7 6 = 1/2, 7 7 = -1/2 • 

Again, according to the extended form of detailed balance, here are two possibilities of directions of 
irreversible reactions 2, 3, 4, 6 and 7 if reactions 1 and 5 are reversible: 2 — >, 3 <—, 4 — 6^,7^ and 

2 3 4 <-, 6 7 

In the third case, let us exclude the coordinate A4 from 72345 ((23)) using vector 7 6 . We get one- 
dimensional vectors 

7 2 = 3, 73 = -3, 74 = 3, 75 = 3 • 

According to the extended form of detailed balance, here are two possibilities of directions of irreversible 
reactions 2, 3, 4, and 5 if reactions 1, 6 and 7 are reversible: 2 — >, 3 , 4 — >, 5 — > and 2 , 3 — k 4 <— , 5 «— . 

In the first and the third cases, there are nontrivial Wegscheider identities for the reaction equilibrium 
constants of reversible reactions. If reactions 1, 2, 3 and 4 are reversible (case 1) then dim A = 2 and the 
basis of A is, for example, A 1 = (2, -1, -1, 0) (2 7i - 72 - 73 = 0) and A 2 = (1, -1, 0, 1) (71 - 72 + 74 = 0). 
The two correspondent Wegscheider identities are: K\ = K 2 K 3 and K1K4 = K 2 (where Ki = k+ /k~^). 

If the reactions 1, 6 and 7 are reversible then dim A = 1 and the basis of A consists of one vector 
A = (1, 1, 1) (71 + 76 + 77 = 0). The correspondent Wegscheider identity is: K\K§K-j = 1. 

If we add one more reversible reaction in cases 1-3 then all the reactions 1-7 should be reversible in 
according to the extended form of detailed balance. 

In this case study, we demonstrated also how it is possible to organize computations and reduce the 
dimension of the computational problems. 

4. Multiscale Degenerated Equilibria 

Let in a system of reversible reactions with detailed balance some k~ — > 0, when the correspondent 
kf remains constant and separated from zero. In this case, some equilibrium activities also tend to zero. 
Indeed, at equilibrium wf = w~ , wj because kj — > 0, hence, wf — > and some of a^ q with a S i > 
also tend to zero due to the generalized mass action law ([5]). Therefore, the irreversible limits of the 
reactions with detailed balance are closely related to the limits when some equilibrium activities tend to 
zero. (For the usual mass action law is sufficient to replace the "activity aj" by the "concentration Cj".) 

In this section we study asymptotics a° q = const x e Si , e — > for various values of non-negative 
exponents 5j > (i = 1, . . . , n). 

There exists a well known way to satisfy the principle of detailed balance: just write k~ = fc+ / K r 
where K r is the equilibrium constant: 

K _ niUKY" 

r nr=i(O a " ' 
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We can define the equ ilibrium constant th r ough the equilibrium thermodynamics as well (see, for example, 



the classical book by iPrigogine fc Defavl (|19 62)). In this case, the principle of detailed balance is also 
satisfied for the mass action law. 

In this approach, we have to group the direct and reverse reactions together. Therefore, m is here the 
number of pairs of reactions, direct + inverse ones. We deal with m + n constants (m rate constants fc+ for 
direct reactions and n equilibrium data for individual reagents: equilibrium concentrations or activities) 
instead of 2 m constants kp. For these m 4 - n co nstants, the principle of detailed balance produces no 
restrictions ( Gorban et all . fl989l : I Yang, et all 2006 ). It holds "by the construction" for any positive values 



of these constants if k~ = kj / K r and the equilibrium constants are calculated in accordance with the 
equilibrium data. 

To transform the conditions of ap — ¥ into irreversibility of some reactions, it is not sufficient to know 
which a° q — > 0. We have to take into account the rates of these convergence to zero for different i. In the 
simple example, A\ ^ A 2 ^ A 3 ^ A l5 if ap 2 — > 0, ai/a 2 — > then in the limit we get the system A\ — > A 2 
(because the Ai/A% equilibrium is shifted to A 2 ), A\ — > A3, A 2 — > A3. For the inverse relations between 
ai and a 2 , 0,2/0,1 — > 0, the limit system is A 2 — > A\ (the A\/A 2 equilibrium is shifted to A\), A\ — > A3, 
A 2 — > A3. For the both limit systems, the equilibrium activities of Ai, A 2 are zero but the directions of 
reaction arc different . 

The limit structure of the reaction mechanism when some of a° q — >• depends on the behavior of the 
ratios a^ q /a^ q . To formalize this dependence, let us introduce a parameter e > and take a° q = const x e 5i . 
At equilibrium, each monomial in the generalized mass action law is proportional to a power of e: 

w c r i+ = k+ const x e^i ; = fc- const x PrA . 

The principle of detailed balance gives: u>° q+ = u>° q ~. Therefore, 

K 



const x e^' 5 ) , (25) 



where 6 is the vector with coordinates Si. 

There are three possibilities for the reversibility of an elementary reaction in asymptotic e — > 0: 

1. If (7r,<5) = then the reaction remains reversible in asymptotic e — ¥ 0. This means that < 
lim(fe]j"/fc~) < 00. Therefore, if one of the reactions persists in the limit then the reverse reaction 
also persists. 

2. If (7,., 6) < then in asymptotic e —¥ can remain only direct reaction. This means that lim(fc~/fc+) = 
0. 

3. If (7 r ,<5) > then in asymptotic e — > can remain only reverse reaction. This means that 
Km(k+/k-) = 0. 

It is possible that (7 r ,<5) = but both fcj, lm± = just because fc+ = and k~ = and not because of 
the equilibrium degeneration. If we delete some irreversible reactions or several pairs of mutually reverse 
reaction then the extended form of detailed balance persists. Therefore, we do not consider these cases 
separately and always discuss the limit reaction mechanisms with max{/cj, lm+ , fcj, lm ~} > 0. 

For each stoichiometric vector 7,. the n-dimensional space of vectors S is split in three sets: hyperplane 
(TV) S) = (reaction remains reversible), hemispace (7^, S) < (only direct reaction remains) and hemispace 
(7r><5) > (only reverse reaction remains). For the reaction mechanism, intersections of these sets for all 
7 r (r = 1, . . . ,m) form a tiling of the n-dimensional space of vectors 5. The intersection of all hyperplancs 
(7 r , 5) = corresponds to the initial reversible reaction mechanism. Other sets from this tiling correspond 
to the reaction mechanisms that are limits of the initial reaction mechanism when some of the reaction rate 
constants tend to zero but the principle of detailed balance is valid. In our study, the exponents Sj should 
be non-negative, hence, we have to study the tiling of the positive orthant 5j > in M. n Description of the 
tiling produced by a system of hyperplanes (7 r , 6) = is a classical problem of combinatorial geometry. 

In the usual linear triangle A\ ^ A 2 ^ A3 ^ A\ we have to consider three hyperplanes in the space 
of exponents S = (61,62,53): 61 = 6 2 ((71, 6) = 0), S 2 = S3 ((72, 6) = 0) and 63 = Si ((73, 6) = 0). At least 
one of the exponents should take zero value to keep the overall concentration in equilibrium neither zero 
nor infinite. Let us take <5i = 0. The hyperplanes turn in the straight lines on the plane (^2,^3): = S 2 
((71, S) = 0), S 2 = 63 ((72,^) = 0) and <5 3 = ((73,(5) = 0). The positive octant on the plane (62,63) is 
split in five sets (A)-(E), that correspond to the limits with some irreversible reactions, and the origin: 

• (A) 6 2 = 0, 63 > 0, Ai ^ A 2 , A 3 -> Ai, A 3 -> A 2 ; 

• (B) 63 > 6 2 > 0, A 3 -> A 2 -> Ai, A 3 -> An 
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Figure 2: Tiling of the positive octant of the plane (62, S3) (5i = 1) that corresponds to seven irreversible limits of the reaction 
mechanism. 



• (C) 63 = S 2 > 0, A 3 ^ A 2 , A 2 -> A u A 3 -> Ar, 

• (D) 82 > 83 > 0, A2 — > A3 — > A\, A2 — ► A\ (this case differs from (B) by the transposition 2 -H> 3); 

• (E) 82 > 0, 83 = A\ ^ As , A2 — > Ax, A 2 — > As (this case differs from (A) by the transposition 
2 «-> 3). 

• The origin corresponds to the fully reversible mechanism. 

For a less trivial example, let us analyze the reaction mechanism from Section 13.41 A\ ^ A2 ^ A3 ^ 
A±, A\ + A2 ^ 2 A3. This is a reversible cycle supplemented by a nonlinear step. 
We join the direct and reverse elementary reactions and, therefore, 



r T 



-1 1-1 
1-1 0-1 
01-12 



(26) 



The columns of this matrix are the stoichiometric vectors j r . 

Let us study the tiling of the positive orthant in R 3 by the planes (7 r , 8) = (r = 1, . . . , 4). First of all, 
it is necessary and sufficient to study this tiling of the positive octants in three planes: 81 = 0, or 82 = 0, 
or 83 = because at least one equilibrium concentration should not tend to zero and, therefore, has zero 
exponent. The symmetry between A\ and A2 allows us to study two planes: 8\ = or S3 = 0. 

On the plane 81 = with coordinates 82, 83 we have four straight lines: (71, 8) = (<5 2 = 0), (72, 8) = 0, 
(82 = 83), (73, 5) — (83 = 0) and (74, 8) = (82 = 28s)- These lines divide the positive octant (82,3 > 0) 
into seven parts (Fig. [2]) and the origin: 

1. (A) 82 = 0,8 3 > 0, A 1 ^ A 2 , As ->■ Ax, A 3 -> A 2 , 2A 3 A 1 + A 2 ; 

2. (B) 8 2 > 0, 83 > 82, A 2 -> A t , A3 -> Ai, A 3 -> A 2 , 2A 3 -> A x + A 2 ; 

3. (C) <5 2 = 83 > 0, A 2 -> Ai, A 3 -> Ai, A 3 ^ A 2 , 2A 3 -> A t + A 2 ; 

4. (D) < 83 < 82 < 283, A 2 -> Ai, A 3 ->• Ai, A 2 -> A 3 , 2A 3 -> Aj + A 2 ; 

5. (E) < 82 = 28 3 , A 2 -> A u As -> Ax, A 2 A 3 , 2A 3 ^ A x + A 2 ; 

6. (F) <5 2 > 2J 3 > 0, A 2 -> A 1; A 3 ->• A x , A 2 -> A 3 , Aj + A 2 -> 2A 3 ; 

7. (G) <5 3 = 0, 82 > 0, A 2 -> A x , A a ^ A 3 , A 2 A 3 , A x + A 2 -y 2A 3 ; 

8. The origin corresponds to the fully reversible mechanism. 

The same picture gives us the plane <5 2 = with coordinates 8x, 83: we need just to transpose the 
indexes, 1h2. 

On the plane 83 = with coordinates 81 , 82 the positive octant is divided into five parts and the origin: 

1. 81 — 0, 8 2 > 0, A 2 -> A x , At ^ A 3 , A 2 ->■ A 3 , Ax + A 2 ->• 2A 3 (this is exactly the case (G) from 
Fig.0; 
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2. < Si < S 2 , A 2 -> Ai, Ai -> A 3 , A 2 -> A 3 , A x + A 2 -> 2A 3 ; 

3. < 6i = S 2 , Ai =± A 2 , A x -4 A 3 , A 2 -> A 3 , Ai + A 2 -> 2A 3 ; 

4. <5i > <5 2 > 0, Ai -s- A 2 , A x A 3 , A 2 ->■ A 3 , A 1 + A 2 -> 2A 3 ; 

5. <5 2 = 0, 5i > 0, Ai A 2 , Ai A 3 , A 2 =± A 3 , A 1 + A 2 -> 2A 3 ; 

6. The origin corresponds to the fully reversible mechanism. 

This approach is equivalent to the previous definition of the extended form of detailed balance based 
on the pathway analysis. Indeed, if the reaction mechanism with some irreversible reactions is a limit of 
the reversible mechanism with detailed balance then it satisfies the conditions of the extended form of 
detailed balance. (This is the direct statement of Theorem Q] proved in Section [3~2| ) To prove the converse 
statement, we have to take a system that satisfies the extended form of detailed balance and to find such 
a set of exponents Si > (i = 1, . . . , n) that the system appears in the limit of a reversible system with 
detailed balance when e — > and a eq = const x e Si . 

Let a system with some irreversible reactions satisfy the extended form of detailed balance. We fol- 
low the notations of Theorem [31 jj (j = l,...,r) are the stoichiometric vectors of the reversible re- 
actions and vi,...,i/ s are the stoichiometric vectors of the irreversible reactions. The linear subspace 
S = span{7i, . . . , 7 r } C K™ consists of all linear combinations of the stoichiometric vectors of the reversible 
reactions. We use notation T7j for the images of Vj in W / S. 

Let kf > (J = 1, . . . , r) be the reaction rate constants for the reversible reactions and qj = > 
(j = 1, . . . , s) be the reaction rate constants for the irreversible reactions. We extend the system by adding 
the reverse reactions with the constants qj > 0. If the extended system satisfies the principle of detailed 
balance then 

k + n a + " 

-^HKT" and ^=11(0^, (27) 

K j i=i q i i=i 



where a t is a point of detailed balance. 

Theorem 4. Let the system satisfy the extended form of detailed balance. Then there exists a vector of 
nonnegative exponents S — (Si) (i = 1, ... ,n) and the family of extended systems with equilibria a° q = a*e Si 
such that condition hold, k^ (j = 1, . . . , r) and qj = q^ (j = 1, . . . , s) do not depend on e, and qj — > 
when e — > 0. 

Proof. If the system satisfies the extended form of detailed balance then the reversible part satisfies the 
principle of detailed balance and, hence, there exists a positive point of detailed balance for the reversible 
part of the system (Theorem [3]): a* > and 

n n 

fc/n«) ari = fc /iK a *) fe - 

i=l i=l 

Let us take a° q = a*e 5i . Due to ([27)) . k^ /kj = const x e^ 5 ) . To keep the kf independent of e, we have 
to provide ("fj,S) = 0. Analogously, q^/qj = const x e^ j ' S K The rate constant q^ should not depend 
on e and qj — > when e — > 0. Therefore, (vj^S) < 0. We came to the system of linear equations and 
inequalities with respect to exponents 

( 7j ,5) = 0(j = l,...,r), {v j ,S)<0(j = l,...,s). (28) 

The solvability of this system is equivalent to the condition 1 of Theorem [3] (see Remark [7]) . To prove the 
existence of nonnegative exponents Si > 0, we have to use existence of positive conservation law: hi > 0, 
(7j, b) = 0, (vj, b) = 0. For every solution 8 of (|2"5)) and any number d, the vector S + db is also a solution 
of (f2"5)). Therefore, the nonnegative solution exists. We proved the theorem and the converse statement of 
Theorem [2 □ 

Proposition 4. Let a system with the stoichiometric vectors 7 S and the extended detailed balance be 
obtained from the reversible systems with detailed balance in the limit a eq = const X — > 0. For this 
system, the linear function (S, c) of the concentrations c monotonically decreases in time due to the kinetic 
equations = Y, s w sls- 



Proof. Indeed, = ~^2 s w s (^ s ,5) (compare to Remark [5)). For the reversible reactions, the sign of 

is indefinite but (7 S , 8) = 0. For the irreversible reactions, we always can take w s = > just by 
the selection of notations. In this case, only fc+ survived in the limit e — > 0, this means that (7^, 8) < 0. 
Therefore, ^) < and it is zero if and only if all the reaction rates of the irreversible reactions vanish. □ 
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So, the vector of exponents 8 defines the (partially) irreversible limit of the reaction mechanism and, 
at the same time, gives the explicit construction of the special Lyapunov function for the kinetic equations 
of the limit system. 

In this Section, we developed the approach to the systems with some irreversible reactions based on 
multiscale degeneration of equilibria, when some ai — > as e Si . We proved in Theorem[4]that this approach 
is equivalent to the extended form of detailed balance based on the pathways analysis or on the limits of 
the systems with detailed balance when some of the reaction rate constants tend to zero. 

5. Conclusion 

The classical principle of detailed balance operates with mechanisms, which consist of fully reversible 
elementary processes (reactions). If such mechanisms have cycles of reactions, each cycle is characterized 
by one Wegscheider relationship © between its rate constants. The number of functionally independent 
relationships is equal to the number of linearly independent cycles, linear or nonlinear. 

In difference from this classical case, we analyzed mechanisms, which may include irreversible reactions 
as well. For such mechanisms we proved an extended form of detailed balance considering the irreversible 
reactions as limits of reversible steps, when the rate constants of the corresponding reverse reactions 
approach zero. The novelty of this form is that the extended detailed balance now is presented as a 
necessary combination of two constituents: 

• Structural conditions in accordance to which the irreversible reactions cannot be included in oriented 
cyclic pathways. 

• Algebraic conditions which are written for the "reversible part" of the complex mechanism taken 
separately, without irreversible reactions, using the classical Wegscheider relationships. 

The computational tools combine linear algebra (some standard tools for chemical kinetics) with meth- 
ods of linear programming. The most expensive computational problem appears when we check the struc- 
tural condition of the extended form of detailed balance. 

Let n be the number of components, and let K™ be the composition space. We consider a system with r 
reversible and s irreversible reactions. Let us use 71, . . . ,j r for the stoichiometric vectors of the reversible 
reactions, f\, . . . , v s for the stoichiometric vectors of the irreversible reactions and Vj for the images of Vj 
in the quotient space M. n /S, where S is spanned by the stoichiometric vectors of all reversible reaction, 
S = span{7i, . . . , 7,.} C R™. The reaction mechanism satisfies the structural condition of the extended 
form of detailed balance if and only if 

^ conv{Fi, . . . ,v s ] . 

We have to check whether the origin belongs to the convex hull of the vectors V\ , . . . , V s . In practice, 
we can always assume that these vectors have exactly known rational (or even integer) coordinates. 

We combined three approaches to study the restrictions implied by the principle of detailed balance in 
the systems with some irreversible reactions: 

1. Analysis of limits of the systems with all reversible reactions and detailed balance when some of the 
reaction rate constants tend to zero. 

2. Analysis of the Wegscheider identities for elementary pathways when some of the reaction rate con- 
stants turn into zero. 

3. Analysis of limits of the systems when some equilibrium concentrations (or, more general, activities) 
tend to zero. 

We proved that these three approaches are equivalent if we take into account not only which equilibrium 
concentrations tend to zero, but the speed of this tendency as well. The various partially or fully irreversible 
limits of the reaction mechanisms are, in this sense, multiscale asymptotics of the reaction networks when 
some equilibrium concentration tend to zero with different speed. 
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